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Abstract. We present high-resolution maps of the CO(l-O) and CO(2-l) emission from the LINER 2 galaxy 
NGC 1961. This galaxy is unusual among late-type (Sc) disk galaxies in having a very large radial extent and 
inferred dynamical mass. We propose a head-on collision scenario to explain the perturbed morphology of this 
galaxy- both the off-centered rings and the inflated radius. This scenario is supported by the detection of a steep 
velocity gradient in the CO (1-0) map at the position of a southwest peak in radio continuum and near-infrared 
emission. This peak would represent the remnant of the disrupting companion. We use numerical models to 
demonstrate the plausibility of the scenario. While ram pressure stripping could in principle be important for 
shocking the atomic gas and produce the striking head-tail morphology, the non detection of this small galaxy 
group in X-ray emission suggests that any hot intragroup medium has too low a density. A prediction of the 
collision model is the propagation of ring waves from the center to the outer parts, superposed on a probable 
pre-existing m — 2 barred spiral feature, accounting for the observed complex structure of rings and spokes. 
This lopsided wave accounts for the sharp boundary observed in the atomic gas on the southern side. Through 
dynamical friction, the collision finishes quickly in a minor merger, the best fit being for a companion with a mass 
ratio 1:4. We argue that NGC 1961 has a strongly warped disk, which gives the false impression of a nearly face-on 
system; the main disk is actually more edge-on, and this error in the true inclination has led to the surprisingly 
high dynamical mass for a morphologically late-type galaxy. In addition, the outward propagating ring artificially 
enlarges the disk. The collision de-stabilizes the inner disk and can provide gas infiow to the active nucleus. 

Key words. Galaxies: individual: NGC 1961 — Galaxies: interactions — Galaxies: ISM — Galaxies: kinematics 
and dynamics — Galaxies: nuclei — Radio lines: galaxies 
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1. Introduction 

The tight correlation between black hole mass and central 
velocity dispersion (Af, — a relation, Merritt & Ferrarese 
120011 Tremaine et al. I2002p offers strong evidence that 
galaxies' stellar spheroids and central black holes are built 



2 



F. Combes et al.: Molecular Gas in NUclei of GAlaxies (NUGA) 



up by mechanisms that are closely linked, if not identi- 
cal. In disk galaxies, black holes of increasing mass can 
scatter more disk stars into a bulge (Norman et al. I1996P 
and drive secular evolution along the Hubble sequence. A 
more generally applicable scenario (holding for both early 
and late-type galaxies), however, centers on the balance 
of nuclear fuelling and feedback. If the relative fractions 
of inflowing gas consumed by nuclear star formation and 
accretion onto a black hole are roughly constant, then the 
M, — a relation is naturally explained. 

Since nearly half of all nearby systems host low lumi- 
nosity active galactic nuclei (AGN: Ho et al. I1997ap . the 
supply of fuel to the Schwarzschild radius (not to mention 
circumnuclear star-forming regions) is still very much an 
ongoing process. Because black holes appear to occupy the 
centers of nearly all galaxies with massive bulges, high- 
resolution studies of the gas dynamics in virtually any 
system with a massive spheroid should therefore shed light 
on the role of fuelling in driving the coevolution of black 
holes and host galaxies. This prospect is the motivation for 
our observations of the Nuclei of Galaxies (NUGA) sam- 
ple at the IRAM Plateau de Bure Interferometer (PdBI: 
Guilloteau et al. [1992]). The NUGA survey is designed 
to explore the molecular gas content and mechanisms for 
nuclear fuelling at high angular resolution; initial observa- 
tions of a sample of twelve nearby AGN (Garcia-Burillo et 
al. I2003P will ultimately be combined with a larger sam- 
ple of both active and quiescent nuclei for the purpose of 
statistical comparisons. In this paper, the twelfth in an 
initial series focused on individual objects, we present the 
PdBI data for the NUGA target NGC 1961. 

Judged by its nuclear properties, NGC 1961 is a rela- 
tively unremarkable low-luminosity AGN. Its optical spec- 
trum is classified as a Low-Ionization Nuclear Emission- 
line Region (LINER) by Ho et al. (|1997b[) : it shows no 
broad-line emission, although its narrow-line profiles do 
show asymmetric blue wings (Ho et al. I1997c"| . NGC 1961 
has attracted significantly more attention for being one 
of the most massive known disk galaxies (Rubin et al. 
I1979p . and moreover unusually massive for its late (Sc) 
Hubble type. This peculiarity cannot be attributed to a 
distance error, since NGC 1961 is associated with a group 
of ten galaxies that all lie within 514kms~^ of each other 
at small separations on the sky (Gottesman et al. I2002p . 
NGC 1961 appears to be impinging on this group from 
the northwest: its large-scale optical morphology shows a 
pronounced linear feature to the southeast which connects 
to an anomalous asymmetric arm winding from south to 
west (Figured]). 

CO emission has already been detected in NGC 1961 
(Young et al. I1995| Komugi et al. I2008p . corresponding to 
a molecular gas mass of 2.7 x VS^^ Mq. The HI content 
of NGC 1961 is substantial (nearly 5 x 10^" Mq of atomic 
hydrogen: Haan et al. I2008p . which suggests a very late 
type galaxy. However, the galaxy has a large radial ex- 
tent, both optically and in HI, apparently much larger 
than expected for its morphological type. The HI mor- 
phology revealed by recent VLA observations (Haan et al. 
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Fig. 1. Blue DSS image of NGC 1961. Note the sharp hn- 
ear feature to the southeast (oriented at 45°), and the 
anomalous arm winding from south to west. The "arm 
features" are tightly wound, and it is not easy to identify 
their winding sense. They might even wind opposite to the 
global sense, at some radii. In any case, the morphology 
is better described by rings. 



120081) is striking because of its head-tail shape, with a very 
sharp boundary in the south. This feature suggests a ram- 
pressure shock, although no diffuse X-ray emission has 
been detected by ROSAT'm. this galaxy group (Mulchaey 
et al. 12003)) . Shostak et al. p982]) first noted the sharp 
H I cut-off displaced with respect to the stellar boundary 
and suggested a ram pressure shock was responsible. This 
linear feature of higher H I column density coincides with 
a higher radio continuum surface brightness at frequencies 
from 1.5 to 15.4 GHz (Lisenfeld et al. [T998)) . 

The galaxy as a whole has an unusually high radio 
continuum luminosity relative to its optical (blue) and far- 
infrared (40 — 120/im) luminosities. Lisenfeld et al. (|1998p 
suggest that a second nucleus in their radio continuum 
map may be the remnant of an infalling dwarf galaxy; 
they attribute the high radio luminosity to cosmic rays 
that have been produced in the putative minor merger. 
The observed morphology can constrain its orbit type: 
the lack of extended tidal tail features in the HI map 
suggests that any recent merger at the very least cannot 
have had a prograde encounter geometry with respect to 
NGC1961's sense of rotation. We will use the molecular 
line observations and the numerical simulations reported 
in this paper to help identify which of these two mech- 
anisms is more plausibly responsible for the large-scale 
asymmetries in NGC 1961. 

We describe our observations in Sj2] and their results in 
fJS] focusing on the m — \ asymmetry and off-centered ring 
morphology evident in the CO maps. A dynamical model 
is proposed and arguments developed in §4] to explain the 
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Table 1. Log of PdBI observations of NGC 1961. For each 
dataset, we list the configuration, the number of antennas, 
and the 3 mm/ 1 mm flux densities adopted for the gain 
calibrator 0716+714 (in mJy). 



Date 


Config 


Afant 


0716+714 


18 Dec 2000 


C 


4 


1.3 / 0.7 


21 Dec 2000 


C 


5 


1.0 / 0.9 


26 Mar 2001 


D 


5 


0.7 / 0.6 


31 Mar 2001 


D 


4 


0.5 / 0.4 


09 Jan 2002 


B 


6 


1.0 / 2.1 


03 Mar 2003 


A 


6 


2.0 / 1.8 



perturbed morphology and asymmetries observed on large 
scales as the result of a recent head-on collision. fj5] de- 
scribes in detail the adopted model, the technique of the 
simulations, and their results. Our discussion of the col- 
lision scenario and its implications for AGN fuelling are 
presented in fJ51 and our conclusions in |J71 Throughout 
this paper, we will use the distance D — 52.4 Mpc spec- 
ified by TuUy (|1988p . and a C0-to-H2 conversion factor 
X = NaJIco{i-o) = 2.2 x lO^" cm'^ (Kkms-i)-i as 
recommended by Solomon & Barrett p99ip . 

2. Observations 

2.1. CO interferometry 

We observed the CO(l-O) and C0(2-l) lines in NGC 1961 
with the IRAM PdBI between December 2000 and 
January 2002. During this period, the array included from 
four to six 15 m diameter antennas; deployment in the 
B, C, and D configurations yielded a total of 35 dis- 
tinct baselines ranging from 24 m to 280 m in length. 
In March 2003, the A configuration added longer base- 
lines, up to 420m for the C0(2-l) line, and provided sub- 
arcsecond resolution ( see documentation on configura- 
tions at http://www.iram.fr). Each antenna was equipped 
with both single-sideband 3 mm and double-sideband 
1 mm SIS receivers that were used simultaneously; receiver 
temperatures were 60-80 K, and corrected outside the at- 
mosphere yielded typical system temperatures of ~ 120 K 
and ^ 400 K at 3 mm and 1 mm, respectively. We deployed 
four correlator modules- giving a total of 580 MHz of con- 
tinuous bandwidth at 2.5 MHz resolution- across each of 
the two fines. The CO(l-O) and C0(2-l) fines were ob- 
served in the upper and lower sidebands relative to the 
3 mm and 1 mm reference frequencies. As the phase (point- 
ing) center for the observations, we adopted a galaxy's op- 
tical position, slightly different from that listed in NED- 
aj20oo = 05:42:04.8 and (5j2ooo = +69:22:43.30 

We calibrated our data using the CLIC routines 
in the IRAM GILD AS software package (Guilloteau & 
Lucas l2000|) . For passband calibration, we used the bright 
quasars 3C111, 3C279, and 3C454.3. Observations of the 
blazar 0716+714 interleaved with our source data every 

^ These coordinates were originally derived relative to the 
position of SN 1998eb. 



~ 30 minutes were used to remove phase drifts and ampli- 
tude gain variations. The fiux density of 0716+714 was in 
turn determined at each epoch by comparisons with the 
calibration sources CRL618 and MWC349; we estimate 
the uncertainties in the flux scales to be ~ 20%. Table [T] 
lists the details of our observations and calibrations, in- 
cluding the fiux densities adopted for 0716+714 at both 
3 mm and 1 mm. The fact that 0716+714 is a blazar ex- 
plains its strong variability. 

Before constructing uv tables from our calibrated data, 
we smoothed them to a frequency resolution of 5 MHz (ve- 
locity resolutions 13.2 km s~^ for the CO(l-O) line and 
6.6kms~^ for the C0(2-l) line). The primary beam sizes 
at the two frequencies are 43 and 21.5"respectively. The 
full datasets include the (on-source, six-telescope) equiv- 
alents of 16.8 hours of 3 mm data and 13.9 hours of 1 mm 
data. In the 3 mm data cube, CO(l-O) line emission is 
detected out to 330 km s^^ blueward and 369kms~^ red- 
ward of the systemic zco — 0.01331 ~ Wgys = 3937 km s~^ 
LSR (see §3.2p . We constructed a 3mm continuum chan- 
nel from the upper-sideband channels outside this velocity 
range, and upon mapping (with natural weighting) discov- 
ered the presence of a 1.3+0.2 mJy continuum source offset 
^ 3'.'4 to the northwest of the nucleus (see i )3.4p . In order 
to avoid contamination of the CO(l-O) line maps, we used 
the GILDAS task UV_SUBTRACT to remove a uv model 
for this continuum source from all of the 3 mm uv data. 
Our maps of the continuum-subtracted line cube used nat- 
ural weighting, which (as for the 3 mm continuum data) 
yielded a synthesized beam of 2'.'73 x 2'.'41 at position angle 
36°. The rms per 13.2kms~^ channel is 1.3 mJy beam~^. 

At 1 mm, no continuum emission was evident (with 
an rms of 0.4 mJy) in a map made from the upper 
(line-free) sideband visibility data, leading us to conclude 
that no subtraction of continuum from our 1 mm line 
data would be necessary. CO (2-1) emission is seen over 
a velocity range of ±257kms~"'^ relative to Vgys, smaller 
than the full CO(l-O) range because gas moving at the 
highest velocities in the CO(l-O) data cube lies outside 
the smaller C0(2-l) primary beam. We used uniform 
weighting to reach an angular resolution of 0'.'91 x 0.'.'84 
(~ 231 pc X 214pc in projection) at position angle 24°. 
The rms per 6.6kms~^ channel is 4.3mJybeam~^. 

We measured the integrated intensity ratio 7?.2i/io — 
-^co(2-i)/^co(i-o) a function of position using maps 
made from uv datasets truncated at an identical radius of 
B/X = 1.6 X 10^ for both lines (i.e., slightly larger than the 
minimum uv radius reached for the C0(2-l) data on the 
shortest projected baselines). The loss of short-spacing 
information means that a map made from the truncated 
CO (1-0) data has a smaller synthesized beam than the 
original version made from the full data; the new C0(2- 
1) map, in contrast, is hardly changed at all. We then 
applied corrections for the primary beam response to both 
maps, convolved the C0(2-l) map to the resolution of the 
(w-y-truncated) CO(l-O) map, and divided them. 
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Table 2. Parameters for NGC 1961 



Parameter 


Value 


Reference 


+20 - 


R.A. (J2000) 


05''42™04.6= 


NED 


+ 15 ■ 


Dec. (J2000) 


69°22'42.4" 


NED 


Vhcl 


3934 km s"^ 


NED 


+ 10 ■ 


RC3 Type 


SAB(rs)c 


RC3 


Diameters 


4.6' X 3' 


LEDA 


c-i +5 - 


Inclination 


47° 


LEDA 


Position angle 


85° 


LEDA 


6 

Q - 


Distance 


52.4 Mpc (1"= 254 pc) 


T88 


(U 
,> 


M(HI) 


4.7 X 10" Mq 


H08 


"o 

-5 - 


Lb 


8.9 X 10" L© 


LEDA 




I/FIR (40 - 120 ^im) 


10.7 X 10" Lq 


IRAS 


-10 - 


R.A. (J2000) 


05'' 42^04.7" 


CO-dyn 




Dec. (J2000) 


69°22'42.3" 


CO-dyn 


-15 ■ 



RC3: de Vaucouleurs et al. (fTMlj) 
T88: TuUy (|1988)) 
H08: Haan et al. ipHnS]) 

CO-dyn: Dynamical center derived from the CO (2-1) kinemat- 



2.2. Optical and infrared imaging 

We have extracted optical and near-infrared images of 
NGC 1961 from the Hubble Space Telescope (HST) archive. 
Images were acquired using WFPC2 through the F547M 
and F606W (- V) filters, and using NICMOS through 
the F160W (- H) filter (see Martini et al. fMJ^ Hunt & 
Malkan[2004j). To show the large-scale optical morphology 
of the galaxy, we have taken red and blue optical images 
from the Digitized Palomar Sky Survey (DSS); the blue 
image is shown in Figure [TJ 

For wide-field near-IR information, we have turned to 
the Two Micron All Sky Survey (2MASS: Skrutskie et al. 
nM7| . Jarrett et al. (^005)) have included NGC 1961 in 
their 2MASS "Large Galaxy Atlas" ; we have downloaded 
their reduced images from the IPAC website. A image 
is shown in Section U) 

IRAC images at 3.6 /xm, 4.5 /im and 8.0 ^m were re- 
trieved from the Spitzer archive. We started with the Basic 
Calibrated Data images, and aligned and combined them 
with MOPEX (Makovoz & Marleau l^UgS]) which accounts 
for distortion and rotates to a fiducial coordinate system. 
1'.'20 pixels were imposed for the final images, roughly the 
same as the original IRAC detector. Significant banding 
was present and this was corrected for by interpolation 
over the affected rows before combination with MOPEX. 
The 8.0 /im IRAC filter tends to be dominated by aromatic 
features of polycyclic aromatic hydrocarbons (PAHs) , but 
also contains a contribution from stellar photospheres. 
The IRAC images at shorter wavelengths (3.6, 4.5 /im) 
can be suitably combined and scaled to estimate the stel- 
lar emission. Following Helou et al. (2004), we obtained a 
non-stellar (dust) image by subtracting from the 8.0 /im 
image a scaled linear combination of the 3.6 and 4.5 /im 
images. The IRAC images were originally published by 
Pahre et al. pOOij) . 
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Fig. 2. CO(l-O) zeroth moment map, with coordinates 
defined relative to the PdBI phase center. The 2.73" x 
2.41" synthesized beam is plotted at lower left; contours 
are {0.04, 0.08, 0.12, 0.16, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 
1.5, 2.0, 2.5} X13.2 Jybeam-^kms-i. 



3. Results 

3.1. CO line morphologies 

We derived zeroth moment maps from our CO data cubes 
using the Groningen Image Processing SYstem (GIPSY: 
van der Hulst et al. 11992^ : emission was integrated over 
the observed line widths after clipping all < 3a pixels 
in the individual channels. Figure [2] shows the result for 
the CO (1-0) line. The global molecular gas distribution 
is dominated by a strong central concentration of dimen- 
sions ^ 10" X 6". This structure is encircled to the north 
and east by a winding gaseous arc containing a number 
of knots, the strongest of which lie ~ 15" east and ~ 7" 
northwest of the nucleus. To the southwest, there is an iso- 
lated concentration of gas whose relationship to the rest of 
the disk is not immediately apparent. An important clue 
to its nature comes from Figure [31 which shows CO (1-0) 
contours overlaid on a greyscale image of 8.4 GHz contin- 
uum emission from Lisenfeld et al. (Ji)98j. Although the 
CO (1-0) primary beam does not extend over the entire 
8.4 GHz map, a comparison in the inner 1' is already 
instructive. The anomalous CO(l-O) emission southwest 
of the nucleus coincides almost exactly with the off-center 
8.4 GHz feature that Lisenfeld et al. (|1998p identify as a 
"second nucleus." This agreement adds some weight to the 
authors' suggestion that the off-nuclear peak represents 
the remnant of a smaller galaxy that has merged with 
NGC 1961 (see also g] below). There is a slight offset be- 
tween the two positions, which is certainly expected, since 
the small companion is likely to be strongly perturbed and 
the CO emission itself is lopsided. 
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Fig. 3. CO(l-O) zeroth moment map contom'S (as in 
Figm'e [2]) overlaid on the 8.4 GHz radio continuum map 
from Lisenfeld et al. (1998) in greyscale. Coordinates are 
defined relative to the VLA phase center adopted for the 
8.4 GHz observations, whose ~ 8" synthesized beam is 
plotted at lower left. CO (1-0) emission coincides with 
the two strongest 8.4 GHz peaks, one in the nucleus of 
NGC 1961 and the other offset to the southwest. 

In Figure m we show an overlay of the CO(l-O) con- 
tours of Figure [2] on the (inner region of the) blue DSS 
image of Figure [T] in greyscale. The anomalous CO(l-O) 
feature to the southwest of the nucleus coincides with a 
region of high extinction in the optical image, suggest- 
ing that the molecular gas structure may lie on the near 
side of the galactic disk. This region also corresponds to 
the possible second nucleus. The eastern section of the 
winding arc also aligns very nicely with a region of high 
extinction, very likely a dust lane, although the northern 
section does not. This contrast offers the first of several 
hints that the CO (1-0) arc does not trace a single spiral 
arm. Instead, as will be clear in the C0(2-l) map, the CO 
emission suggests a distorted ring-like structure. 

The reliability of the CO (1-0) line luminosity as an ex- 
tragalactic mass tracer (e.g., Dickman et al. 11986) Solomon 
& Barrett I1991|) allows us to convert the line fluxes of 
various structures in the CO(l-O) map into molecular gas 
masses. Including a factor of 1.36 to account for helium, 
the gas mass bound into molecular clouds can be calcu- 
lated as 

= 3.23 X 10^ f 
Mq VJykms-i/ 

^(2.2 X 1020) ^ (52.4Mpc) 

After correcting Figure [5] for primary beam response, we 
estimate the masses of the central disk, and of the CO 
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Fig. 4. PdBI CO(l-O) zeroth moment map contours (as 
in Figure [2|) overlaid on the blue DSS image in greyscale. 
Coordinates are defined relative to the PdBI phase center. 

ring including the anomalous southwestern emission, to 
be 1.4 X 10^° Mq and 1.1 x 10^° Mq, respectively. The 
molecular mass corresponding to the southwestern emis- 
sion alone is 1.5 x 10^ Mq. Summing the CO(l-O) line flux 
over the galaxy's observed area gives a total molecular 
mass of Mgas = 2.5 x 10^° Mq. This molecular component 
is more than an order of magnitude more massive than in 
any NUGA galaxy studied so far. 

Figure [5] shows the zeroth moment of the CO (2-1) 
line. Due to primary beam attenuation, only the strongest 
(~ 7" northwest) knot in the northern portion of the 
CO(l-O) arc is detected in C0(2-l). However, the cen- 
tral concentration in the CO (1-0) map is clearly resolved 
into a more compact structure encircled to the south by 
a clearly continuous but strongly asymmetric spiral arm 
or arc. This southern arm coincides with an optical dust 
lane, which is better seen on the C0(2-l)/unsharp mask 
NICMOS image overlay shown in Figure [6l The southern 
arc winds up toward the north in a ring, and lacks the 
constant-sign pitch angle that an arm would have. This 
overlay has been done by making the dynamical CO (2-1) 
map coincide with the peak in the NICMOS image (cor- 
responding also to the peak in the C0(2-l) image). The 
NICMOS image in Fig [S] (top) shows some evidence of a 
nuclear stellar bar of diameter roughly 1 kpc (r^4 arcsec) , 
oriented along PA — - 50°. The signature of the nuclear 
bar is clearly seen in the F160W radial brightness pro- 
file (not shown) as an ellipticy peak over a constant PA 
~ —50°. This is clearly a nuclear bar, since its elongation 
and position angle is different from the rest of the disk. 

In Figure[7l which shows C0(2-l) contours overlaid on 
CO(l-O) greyscale, it appears that the tip of the southern 
arm in the C0(2-l) map is more or less continuous with 
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+4 +2 -2 
relative R.A. ["] 

Fig. 5. C0(2~l) zeroth moment map, with coordinates 
defined relative to the PdBI phase center. The 0.91" x 
0.84" synthesized beam is plotted at lower left; contours 
are {0.12, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 1.5, 2.0, 3.0, 4.0, 
5.0} x6.6 Jybeam"^ kms~^. 



the northern section of the CO(l-O) arc. It is obvious, 
however, that the complex structure of arcs does not cor- 
respond to a simple and classical m = 2 spiral structure. 
Instead, the pitch angle of the arcs is very low, and they 
wind up over almost a turn, more closely resembling a ring 
morphology. The CO is distributed in 2 or 3 embedded, 
offset rings, suggesting an m = 1 asymmetry. This ringed 
C0(2-l) + CO(l-O) structure would then continue to be 
embedded in another ring structure seen on large scales in 
the optical. The external ring is much more asymmetrical 
and off-centered. These features provide strong clues as 
to the nature of the event that produced the large-scale 
asymmetries summarized in iJT] Although it is possible 
to explain the small-scale lopsideness by internal mecha- 
nisms only (e.g.. Jog & Combes r2009p . the external ring's 
off-centering suggests a galaxy collision. 

The innermost C0(2-l) contours reveal a double- 
lobed structure of dimensions ~ 3" x 2" that is not aligned 
with the near-infrared nuclear bar in the HST-NICMOS 
image of Figure [6] This structure could correspond to the 
gas response to the nuclear bar, although this is not the 
only possible explanation in all cases (Schinnerer et al. 
mm Baker & ScoviUe T^SS). The peak in the C0(2-l) 
emission occurs in the southern lobe at an offset of (—1", 
— 1") with respect to the PdBI pointing position (Figure 
[5]), and coincides with the NED optical center (see Table 

El). 




-a 




Fig. 6. Top C0(2-l) zeroth moment contours (as in 
Figure [5]) overlaid on HST F160W colourscale. Bottom 

Same contours superposed on the Unsharp Mask version 
of the NICMOS image. 



3.2. CO line kinematics 

We have also used GIPSY to analyze the CO velocity fields 
in NGC1961. Figures [8] and [9] show the first moments of 
the CO(l-O) and 00(2-1) lines, constructed with clipping 
of (respectively) < 5(t and < 3cr pixels in the original ve- 
locity channels. As expected given the orientation of the 
galaxy's isophotal major axis on the sky, the predominant 
velocity gradient in both first moment maps runs in the 
same east-west direction; the ascending node has a posi- 
tion angle ~ 263°. The larger amount of dust seen towards 
the southern part of the disk, both in the DSS image on 
large scales and in the HST image on small scales, estab- 
lishes that the near side is south and the far side is north. 
The receding velocities in the western half of NGC 1961 
then imply that most arcs are trailing, but some are also 
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Fig. 7. C0(2-l) zeroth moment contours (as in Figure[2]) 
overlaid on CO(l-O) zeroth moment greyscale. 




+25 +20 +15 +10 +5 -5 
relative R.A. ["] 

Fig. 8. CO(l-O) first moment map, with coordinates de- 
fined relative to the PdBI phase center. The synthesized 
beam at lower left is as in Figure [21 contours are in steps 
of 20kms~^ relative to Vsys, with negative (positive) ve- 
locities indicated by black (white) curves. 



leading. As noticed in the previous section, many arcs 
are part of large-scale rings, with almost zero pitch an- 
gle. When spokes are seen between two rings, these are 
unambiguously trailing. 

In contrast to the fairly regular velocity field in the 
main disk, however. Figure [8] shows an obvious disconti- 
nuity to the southwest of the nucleus at the location of 
the putative second nucleus. If considered as a distinct 
structure, the velocity gradient of ~ 200kms~^ along 
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Fig. 9. C0(2-l) first moment map, with coordinates de- 
fined relative to the PdBI phase center. The synthesized 
beam at lower left is as in Figure EJ contours are as in 
Figure [H 



a kinematic major axis of ^ 135° would imply an en- 
closed dynamical mass of ~ 3.0 x 10^ Mq within a ra- 
dius of 5" ~ 1.3 kpc, assuming an edge-on system. This 
is 3 times smaller than the central molecular gas mass 
estimated in i i3.1l above. However, the dynamical mass 
could be higher, if a more face-on inclination is assumed. 
Because the range of velocities seen in the southwestern 
peak is reasonably continuous with those seen in the ad- 
joining region of the disk, we cannot completely rule out 
the possibility that it is in the disk plane, but experiencing 
rather strong and irregular streaming motions. In Figure 
|8l two other kinematically perturbed regions can be seen 
at Aa,AS) = (-17", -2'.'5) and Aa, A(5) = (12'.'5,0"), 
with a 100 kms^^ gradient over 3". Both correspond to 
large molecular condensations, that provide their own self- 
gravity, in addition to the large-scale velocity gradient. 
These do not require extra mass. 

At the higher resolution afforded by Figure[9l it is clear 
that the true dynamical center of NGC 1961 lies between 
the two lobes seen in Figure[5l rather than coincident with 
either of them, with the optical peak, or with the coordi- 
nates we adopted as the phase center for the PdBI ob- 
servations. A fit to the C0(2-l) velocity field locates the 
dynamical center ~ 1" south and ~ 0.5" west of the phase 
center, i.e., at a position aj2ooo = 05:42:04.7 and (5j2ooo = 
-1-69:22:42.3. With respect to this center, we have derived 
circular velocity curves for NGC 1961 in both CO lines. 
Using the GIPSY task ROTCUR in automated mode, we 
have derived fits to the CO(1~0) and C0(2-l) velocity 
fields depicted in Figures [5] and ^ the results are pre- 
sented as dashed lines in Figure [TUl The effects of beam 
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ROTATON CURVE OF NGC 1961 



-■■2 1^ II; li Z.i AL L\ J!j AM JU 



Fig. 10. CO-derived rotation curves for NGC 1961. 
CO(l-O) and C0(2-l) curves are plotted in black and 
gray, respectively; fits to the two-dimensional velocity 
fields with ROTCUR and to the three-dimensional data 
cubes with INSPECTOR are plotted as dashed and solid 
lines, respectively. None of the curves are corrected for 
inclination. 



smearing in reducing the dv/dr slope of the CO (1-0) rota- 
tion curve relative to the higher-resolution C0(2-l) data 
are immediately apparent. As an alternative, we have also 
used the interactive GIPSY task INSPECTOR to derive 
rotation curves from the full data cubes. These results are 
shown as solid lines in Figure [TOl Because the fits to the 
full data cubes are less strongly affected by beam smear- 
ing, they yield higher peak velocities at small radii, and 
thus convey a more accurate impression of the high mass 
densities in the galaxy's inner arcseconds. 

The fuU complexity of the velocity field in NGC 1961 
can perhaps best be appreciated from inspection of the 
four panels of Figure fTTI Here we have extracted position- 
velocity cuts along the kinematic major (position angle 
263°) and minor (position angle 173°) axes. The spatial 
widths of the cuts are 50" for the CO(l-O) cube and 
20" for the C0(2-l) cube. Along the minor axis, we see 
substantial central velocity dispersions and indications of 
streaming motions within the inner 3" ^ 763 pc. Along the 
major axis, the CO(l-O) cut reveals that both the west- 
ern section of the arc (at the most positive X-ofFsets) and- 
especially- the eastern concentration (at the most nega- 
tive offsets) are moving at velocities faster than would be 
expected if motions at smaller radii were straightforwardly 
extrapolated. The C0(2-l) major-axis cut is markedly 
asymmetric and turns over towards the east. This pattern 
can be explained either by strong streaming motions in 
the central 6", or by a scenario in which the eastern 
CO (1-0) concentration has its own sense of motion and is 
not following circular rotation in the disk plane. 



(a) PA 263°-- 




(b) PA 173° 



6> 
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Fig. 11. Position-velocity cuts along the kinematic major 
(P. A. - 263°, West is left) and minor (P. A. ~ 173°, 
South is left) axes of the CO(l-O) (a & b) and C0(2-l) (c 
& d) data cubes. Contours are in steps of {3, 6, 12, 24,...} 
xlcr. 



u 




Fig. 12. Map of the CO(2-1)/CO(1-0) ratio obtained af- 
ter smoothing the CO (2-1) map to the resolution of the 
CO(l-O) map. The contours are from 0.8 to 2.0 in steps 
of 0.2. The ratio increases from 0.7 in the center to 1.5 in 
the ring. 
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Fig. 13. 3 mm continumm contours overlaid on the C0(2- 
1) zeroth moment map in grayscale. The resolution of the 
3 mm continuum map is 2.73" x 2.41" at position angle 
36°. Contours are at 0.6 and O.TmJy beam~^, starting at 
3(j (rms is 0.2mJybeam~^). The dynamical center de- 
rived from the C0(2-l) velocity field in g3?2] is at RA= 
05''42™04.7'^ and Dec= 69°22'42.3". 

3.3. CO (2-1) /CO (1-0) intensity ratio 

Figure [12] shows our measurement of the intensity ratio 
7^21/10 for regions where both CO(l-O) and C0(2-l) maps 
have surface brightnesses > 2a. Over most of the disk of 
NGC 1961, we measure 7?.2i/io to be 0.7— 1.0, close to the 
mean ratio 0.89 ±0.06 measured for local galaxies (Braine 
& Combes ll992p and consistent with optically thick emis- 
sion from thermally excited lines. However, contrary to 
most galaxies, where the ratio decreases from the center 
outwards, here the ratio increases from 0.7 in the very cen- 
ter to 1.5 in the C0(2-l) ring of radius 5 "= 1.3 kpc. The 
high excitation in the ring could be due to more intense 
star formation. 

3.4. Millimeter continuum emission 

NGC 1961 is one of a growing number of nearby AGN 
which reveal continuum emission at millimeter wave- 
lengths (other examples include NGC 1068, NGC 3031, 
NGC 3147, NGC 3718 and NGC 4579: Heifer et al. 'WOT, 
Krips et al. 120051 1200 6 ). Our non-detection of 1mm con- 
tinuum emission (with rms of 0.4 mJy) excludes dust as 
a possible explanation, except in the unlikely event that 
the 0.91" X 0.84" beam at 1mm has resolved out a struc- 
ture that was unresolved at the 2.73" x 2.41" resolution 
of our 3 mm map. If the 1.3 mJy source detected at 3mm 
was coming from dust emission, given typical dust tem- 
peratures, a 1mm intensity of at least 5 mJy would be 
expected, which would have been detected at 12(t. A flat 
spectrum can also be ruled out at 3a. We conclude that 




Fig. 14. Short-wavelength dust (non-stellar) emission 
from the Spitzer IRAC images (see Sect. 2.2). North is 
up. East to the left; the field-of-view is 4.2 arcmin on a 
side. 



the continuum emission must be nonthermal in origin, im- 
plying that the source of ionizing photons responsible for 
the galaxy's optical LINER spectrum is more likely to be 
accretion than star formation. Figure [T3] illustrates the 
^ 3.4" offset of the continuum peak to the northeast of 
the dynamical center estimated from our fit to the C0(2- 
1) velocity field in ^3.21 Although the ~ 8" resolution of 
the Lisenfeld et al. (|1998p radio continuum data are not 
sufficient to permit a detailed comparison, we believe this 
discrepancy is real; however, we do not view it as evi- 
dence for the migration of a supermassive black hole away 
from the center of the galaxy's potential. Instead, as is 
also the case for much of the 3 mm continuum emission 
in NGC 1068 (Schinnerer et al. [2000)) . we might be see- 
ing offset jet emission. The implied emergence of the jet 
from the nucleus at an angle 9 < 90° with respect to the 
disk plane would not be exceptional, given the skewed con- 
figurations inferred for many other nearby low-luminosity 
AGN (Kinney et al. l2000p . At very small scales, centimeter 
VLBI observations provide evidence of extended emission, 
which could be the base of the jet (Krips et al. 12007^ . 

4. A possible minor merger 

NGC 1961 is the largest galaxy in a small group that con- 
tains only spirals, with a velocity dispersion of 182 km s^^ 
too low for the group to contain hot and diffuse X-ray 
emitting gas (Pence & Rots I1997p . The absence of an 
intragroup medium (IGM) does not favor a scenario in 
which ram pressure has a large effect on the perturbed 
morphology. On the other hand, there are no close com- 
panions presently near NGC 1961, although there are sev- 
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eral group galaxies at distance larger than 100 kpc, mainly 
to the south-east (Shostak et al. I1982p . The only plausi- 
ble scenario involving a galaxy interaction is thus a minor 
merger, and there are indeed signs of a possible satellite 
remnant within the confines of NGC 1961. 

The morphological and kinematic information we have 
gleaned from our CO(l-O) and C0(2-l) observations of 
NGC 1961 has provided a set of important clues which we 
can use, in conjunction with data at other wavelengths, to 
try to understand the galaxy's recent dynamical history: 

1. Most of the molecular gas within a radius of 20" ^ 
5.1 kpc is distributed in a disk whose velocity field is 
reasonably regular. At a radius of 21.5" ~ 5.5 kpc, 
however, we see an isolated structure with Mgas 
1.5 X 10^ AIq and (if self-gravitating) M^yn ^ 3 x 
10^ Mq at least, or higher according to its inclina- 
tion, roughly coincident with the location of a peak 
in 8.4 GHz continuum. 

2. The morphologies of the molecular gas structures are 
quite asymmetric, with a one-armed spiral winding 
through roughly 270° in azimuth, or more likely (since 
the pitch angle is close to zero) a molecular ring. The 
CO(l-O) reveals a ring of diameter 35"= 8.9 kpc, and 
the CO (2-1) reveals an asymmetrically embedded ring 
of diameter 9"= 2.3 kpc. A hint of a nuclear gas bar is 
visible in the high- resolution CO (2-1) map, with a di- 
ameter of 0.7 kpc, and a PA of 15°. There is a nuclear 
stellar bar of roughly 1 kpc in diameter, revealed by 
HST in the near- infrared, with a PA of -50°. We note 
that the two CO rings at radii of roughly 1 kpc and 4 
kpc could correspond to the ILR and OLR of this bar, 
if it had a pattern speed of ~ 120kms^^ kpc^^ (from 
Fig.ini) . 

3. The southwestern CO (1-0) emission coincides with a 
peak in the near infrared map, at about 5 kpc from 
the center of NGC 1961. This NIR secondary peak is 
split in three smaller components, suggesting the dis- 
ruption of a stellar system. This splitting could also 
be partially due to dust, although not likeky, given the 
moderate average gas column density in the region. 

We can use these facts to help evaluate the merits of the 
minor merger and ram pressure stripping scenarios which 
have been proposed to account for the large-scale asym- 
metries seen in NGC 1961. An initial, obvious point is that 
because these asymmetries are seen in stars as well as in 
gas (e.g., Figure [T]), the classical version of ram pressure 
stripping conceived as acting only upon gas (Gunn & Gott 
I1972|) will not be sufficient. A more nuanced stripping sce- 
nario would resemble the recent simulations of VoUmer et 
al. (200T) and VoUmer ([20031) . who have shown that it 
is possible for a galaxy's stellar potential to respond to 
a strong perturbation induced by the gas ram pressure. 
However, in the case of NGC 1961, there appears to be 
too little intragroup medium to allow for this. If there 
were sufficient hot IGM, the geometry of the H I tail would 
correspond nicely to what is expected from ram pressure. 
The buildup of HI and enhancement of radio continuum 



emission in NGC 1961's southeastern ridge could then be 
accounted for by an edge-on stripping geometry. For such 
a geometry, simulations predict that a galaxy encounter- 
ing an intragroup medium for the first time should have its 
H I extended on the side where the galactic rotation adds 
to the velocity relative to the IGM and the gas is more 
stripped on the other side (see the explanation of this 
paradox in Roediger & Bruggen |2006l Jachym et al. 120071 
[20091) . For NGC 1961, which is falling into its group from 
the northwest, we see an extended HI structure on the 
northern side of the galaxy, coming from the west (Haan 
et al. I2008P . For ram pressure to produce this configura- 
tion, the velocity of NGC 1961 relative to its group should 
be blueshifted towards us (i.e. the galaxy should be located 
behind the main group). Thus the ram-pressure scenario 
is plausible from a geometrical point of view, the only 
problem is the insufficient IGM. 

In the absence of IGM, we arrive therefore at the con- 
clusion that a tidal encounter or minor merger is the most 
likely explanation for NGC 1961 's current state of disrup- 
tion. The galaxy interaction cannot be a prograde grazing 
passage, as in this case extended tidal HI tails will form, 
which are not observed. The interaction must be either 
retrograde, or nearly head-on. The rings we have detected 
in the CO maps favor a head-on encounter, with a small 
impact parameter. 

Since there is no obvious interaction partner remain- 
ing on the scene, it is also reasonable to conclude that 
the smaller galaxy has already merged with NGC 1961; 
this hypothesis has the advantage of naturally account- 
ing for the distinct kinematic axis of the southwestern 
peak in CO (1-0). We can expect a priori that the merger 
mass ratio should be of order 1:5, to account for the sub- 
stantial (but not total) disruption of the original disk of 
the smaller galaxy. These considerations together provide 
us with starting conditions to develop a detailed minor 
merger model. 

The Spitzer IRAC images bring new insight on the true 
morphology of NGC 1961. As can be seen in Figure[T4l the 
dust emission clearly reveals two off-centered rings, in a 
morphology similar to that of the Cartwheel galaxy (e.g. 
Appleton & Struck-Marcell fT996p . In particular, the ring 
morphology is closely related to spokes in the disk: multi- 
ple spokes link the inner ring to the outer off-centered ring, 
starting nearly tangential to the inner ring and ending per- 
pendicular to the outer ring. These specific details are not 
so apparent in the optical image, which is confused by 
star-formation spots and extinction, although the spokes 
still can be seen. The spokes bridging the gap between off- 
centered rings are typical of a head-on collision in a spiral 
galaxy: either the latter possesses a spiral structure with 
a large pitch angle before the collision, and this structure 
is superposed on rings that propagate outwards from the 
center during the collision (e.g.. Block et al. '2006), or the 
spiral structure forms through gas instabilities during the 
collision (e.g., Hernquist & Weil 1 19931) . 

The ring and spokes morphology of even the cen- 
tral stellar and ISM components in NGC 1961 suggests 
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Fig. 15. The Spitzer 8 /im contours are superposed on 
the color map of the HI surface density from Haan et 
al (PPDSD . 
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Fig. 16. 2MASS X-band image of the central 76.6" x 
55"of NGC1961. The cross- hair indicates our adopted 
phase center for the PdBI observations (the SE star is 
at 29.6"from this center, and can serve as a reference for 
the resolution of the image. North is up, and east is to the 
left). The NIR images all reveal a high stellar density at 
about 24"to the southwest, which is split in three pieces. 
This location corresponds to the high velocity gradient de- 
tected in the CO(l-O) map, and suggests the presence of 
a small disrupting satellite. 



a nearly head-on encounter. Of course, this could also be 
superposed on some ram-pressure phenomenon. However, 
the superposition of the H I and IRAC PAH-dust images 
in Figure [12] shows a rather good correspondence, sug- 
gesting a dominant gravitational phenomenon, which will 
similarly drag the different ISM components, diffuse or 
clumpy. The stellar component, traced by the IRAC 3.6/xm 
image, also follows the sharp southern arc (Pahre et al. 

mm . 

The possible existence of the satellite remnant, sug- 
gested already by Lisenfeld et al. (|1998p from the radio 
map, is also supported by the near-infrared images that 
trace the stellar density. In the 2MASS images (see Figure 
there is clearly a second nucleus, which appears split 
into 3 clumps. Its average distance from the center of 
NGC1961 is 24 "or 6.1 kpc. This feature cannot be seen 
in the HST NICMOS images, since it falls outside of the 
field of view. 

The NOT image of NGC 196l[l constructed by J. 
Naranen & K. Torstensson in 2004 by combining B, V, 
and R exposures shows that the satellite nucleus is par- 
tially hidden by a dust lane, and is split into 2-3 clumps. 
It is clear in this image that star formation is strongly en- 
hanced in the rings, which are delineated by multiple blue 
hot spots. 

Summarizing, the evidence that NGC 1961 has expe- 
rienced a head-on collision followed by a minor merger 
includes: 



^ see the web site http://www.not.iac.es/general/photos/ 
tronomical / extragalactic / 



1. the off-centered ring- like morphology (recalling that 
the main arcs/arms have no obvious sense of wind- 
ing), observed in many components and wavelengths, 
including in CO and in dust emission (e.g.. Fig. [H|) : 

2. the very sharp southern edge of the disk observed in 
the gas and dust component, in particular in the HI 
21cm map (see Figure [15]), which could be the tracer 
of an expanding ring wave; 

3. the existence of "spokes" linking the ring features, with 
a large pitch angle, particularly visible in the Spitzer 
dust image (Figure [T4|: and 

4. the presence of a perturbation in the CO velocity field 
5 kpc from the center, coinciding with a secondary stel- 
lar peak and a radio continuum peak. 

5. Simulation of the head-on collision 

To test the plausibility of a head-on encounter leading to a 
minor merger in NGC 1961, we performed a self-consistent 
simulation of the collision. 

5.1. Galaxy model and numerical techniques 

The simulations are 3D TV-body with stars and gas but not 
including star formation. They are fully self-consistent, 
with a live dark halo, and use a Particle-Mesh code based 
on FFT with a useful grid of 128'^ (the algorithm of James 
ri977 allows us to suppress the Fourier images). The grid 
size is (72kpc)'^, and the model galaxy is truncated at an 
initial radius of 22 kpc. The softening is equal to the size 
of a cell, i.e., 560 pc. The time step is IMyr. The gas is 
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represented by sticky particles, and a total of 2.4 10^ parti- 
cles are used. The non-dissipative components, stars and 
dark matter, are represented by 2 million particles, and 
the gas component by 4 10^. First tests for the simula- 
tions were done with 10 times less particles, and the total 
number were used only when a reasonable fit was found. 
The comparison between simulations with different par- 
ticle numbers reveal that the main features (propagating 
ring waves, off-centering) are a fundamental and persist- 
ing result, while secondary ones, like the exact positions 
of spokes, or the contrast of the arms/rings, are changing. 

The stellar component is composed of a bulge and 
a disk. The bulge is initially distributed as a Plummer 
sphere, with a potential: 



Mr) = - 



GMh 



(2) 



where Mf, and rt, are the mass and characteristic radius of 
the bulge, respectively (see Table [S]). 

The stellar disk is initially a Kuzmin-Toomre disk of 
surface density 



E(r) = S„(l + rVr^)-^/^ 



(3) 



with a mass M^, and characteristic radius r^. It is initially 
quite cold, with a Toomre Q parameter of 1. The dark 
matter halo is also a Plummer sphere, with mass Mdm 
and characteristic radius tdm- The initial conditions of 
the run described here are given in Table [31 

Table 3. Initial conditions 



Component 


radius 


Mass 


Mass fraction 




[kpc] 


[Mq] 


[%] 


Bulge 


0.5 


4.el0 


13.3 


Disk 


6.5 


g.elO 


30. 


Halo* 


15. 


15.el0 


50. 


Gas 


8. 


2.el0 


6.6 


Companion 


2.5 


7.el0 


23.3 



* Halo mass inside 22 kpc radius 

The mass fractions are normalized to the total mass of 
NGC 1961, 3 X 10" Mq inside 22 kpc 



The gas is treated as a self-gravitating component 
in the iV-body simulation, and its dissipation is repre- 
sented by a sticky particle code, following Combes & Gerin 
(1985). We varied the initial gas-to-total mass ratio be- 
tween 4 and 8% to control the stability of the disk prior 
to the collision. The mass of one gas particle therefore 
varied between 3 x 10** and 6 x lO'*Af0. This is meant 
to represent essentially the molecular gas. The initial gas 
distribution is also a Kuzmin-Toomre disk, but truncated 
at 34 kpc, and with a characteristic radial scale of 8 kpc. 
Initially, its velocity dispersion corresponds to a Toomre 
Q parameter of 1. The gas clouds are subject to inelastic 
collisions, with a collision cell size of 560 pc (the volume 
within which particles are selected to possibly collide). 



This corresponds to a lower limit for the average mean 
free path of clouds between two collisions. The collisions 
are considered every 10 Myr. In a collision, the sign of the 
relative cloud velocities is reversed and the absolute values 
are reduced: relative velocities after the collision are only 
fei times their original value, with the elasticity factor 
fei in the radial direction controlling the dissipation rate. 
The fei parameter must be adapted to the particle num- 
ber, so that a given gas particle experiences a reasonable 
loss of energy in a rotation. This is calibrated in simula- 
tions of the isolated galaxy, where there is an equilibrium 
between the dissipative cooling, and the dynamical heat- 
ing due to spiral arm formation. The fei parameter has 
been fixed to 0.65 for low particle numbers, and 0.85 for 
high particle numbers (10 times more particles). The anal- 
ogous coefficient in the tangential direction is fixed to 1, 
to conserve angular momentum in each individual cloud 
collision. All gas particles have the same mass (5 x 10** Mq 
in the standard model) . The particle mass for the dissipa- 
tionless components, stars and dark matter, is identically 
1.5 X IO^Mq. 

The rotation curve corresponding to the standard run 
is plotted in comparison to the data points in Fig. [TT] All 
other runs have similar rotation curves. The parameters 
varied were essentially the geometry of the interaction, the 
mass of the companion, and the stability of the disk before 
the interaction. 

One of the most surprising characteristics of NGC 1961 
is its high dynamical mass relative to its late type. This 
problem is alleviated by the choice of a higher inclination 
of the main disk with respect to the plane of the sky. It will 
be shown in the next section that the outer disk is warped, 
making the galaxy artificially appear more face-on. Here 
we adopt a high inclination of 65-70° for the initial disk 
orientation, and hence for the calculated rotation curve. 
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Fig. 17. Rotation curve and derived frequencies $7,57 — 
k/2 and f2 + k/2 adopted for the simulations presented 
here. The triangles are the CO rotational velocities; the 
squares are from the HI data (Haan et al 2008), with an 
adopted inclination of the galaxy main plane of 65°. 
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The companion is represented by a Plummer poten- 
tial, where radius and mass are given in Table [3l The 
collision corresponds to a minor merger with a mass ratio 
around 1:4. Since we want to explain the morphology of 
NGC 1961, and have no constraint on its satellite, the com- 
panion is simulated as a rigid body. Its orbit is computed 
from the forces exerted on the satellite by all the particles 
in the main galaxy, so that most of the dynamical friction 
is taken into account. The contribution to dynamical fri( - 
tion from the companion disruption is not considered, and 
could decrease the merging time (see Prugniel & Combes 
I1992p . Given the many free parameters, we do not aim to 
reproduce exactly the whole collision or the exact position 
of the companion now, but just to demonstrate that this 
type of collision is able to provide the right order of mag- 
nitude for the perturbation, and morphologies similar to 
those observed. 

5.2. Simulation results 

About a dozen runs have been carried out, to test the 
various geometrical parameters of the encounter and also 
to select the initial morphology of the galaxy. Indeed, a 
bar and a spiral wave develop in the disk during the first 
period when the galaxy simulation is run in isolation, and 
the starting time of the collision plays an important role 
in determining the final morphology. Similar experiments 
were carried out to simulate the head-on collision of the 
Andromeda galaxy: a small satellite (possibly M32) could 
have triggered the propagation of two ring waves in the 
disk, which superposed on the pre-existing spiral structure 
give the complex and perturbed morphology of M31 (cf 
Block et al [2006)1 . 

The orbit of the companion is defined by its energy and 
its impact parameter. The energy is computed from the 
equivalent two-body problem, assuming that all the mass 
of NGC 1961 is concentrated in a point (of mass M). The 
energy per unit mass of the reduced particle is then: 

E ^^V^ ~GM{1 + fi)/R (4) 

where /i is the mass ratio of the satellite to the main 
galaxy, and R and V are the relative distance and velocity 
between the two. 

The initial distance of the companion is 70 kpc, with 
a dominant component perpendicular to the plane of the 
main galaxy. The kinetic energy and relative velocity are 
then obtained, assuming either a "hyperbolic" type of or- 
bit with positive energy, or an "elliptic" type with nega- 
tive energy. The dynamical friction then reduces the rel- 
ative energy from its initial value, and all orbits end up 
as bound. The best run (presented here) has initially a 
negative energy, of —0.87 x lO'' km^ s~^. 

The various experiments revealed that the impact pa- 
rameter must lie between 4 and 10 kpc; it can be no larger 
than 10 kpc if the ring waves characteristic of a head-on 
collision are to appear, but no smaller than 4 kpc if the col- 




Fig. 18. Particle plots from the simulation of the col- 
lision, projected to be similar to the observations (incli- 
nation 70°): Left: the stellar component. Right: the gas 
component. The two snapshots correspond to 320 and 360 
Myr (top and bottom respectively) . The size of each pro- 
jected field is 63 kpc. 



lision is to generate sufficient asymmetry and off-centering 
of the rings. 

During the head-on collision, part of the outer disk is 
raised outside the plane in a huge warp. This global motion 
should be even more visible in the outer gas disk, and we 
suggest that the conspicuous northwestern features seen 
in the H I map, forming a fan with low-projected radial 
velocity along the line of sight, is a consequence of this 
warp. The deformation of the plane to a more face-on 
orientation has thus deceived us as to the true inclination 
of the main plane and the inferred dynamical mass. 

In the run presented in a few snapshots in Fig. 1181 and 
[T51 the initial galaxy model was run in isolation for 280 
Myr, and formed a barred spiral. Then the companion 
was introduced at 70 kpc, and the collision was followed 
during the subsequent 440 Myr. At this epoch, estimated 
to give the best fit to the observations, the ring waves 
have already begun to propagate outwards, and it is easily 
seen how the apparent radius of the galaxy inflates in the 
figures. After 1.5 Gyr of evolution, the galaxy could then 
become a giant LSB galaxy, as proposed by Mapelli et al. 

poosli . 

A density map of the stellar and gas component is 
plotted in Figure [^D] and the corresponding velocity fields 
in Figure [2T] It is possible to see a small accumulation 
of particles accreted by the companion on its way back 
to merge with the main galaxy. Since the companion's 
position is not fit to what is observed, the velocity gradient 
observed around it is located too far north; however, the 
physics of the satellite were not actually addressed in these 
simulations. It is not even known whether the molecular 
gas observed in this region is coming from the main galaxy 
or from the merging satellite. 
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Fig. 19. Same as Fig. [HI but now the snapshots corre- 
spond to 400 and 440 Myr (top and bottom respectively). 
Note how the radius of the galaxy inflates while the ring 
wave propagates outward. 



stars t=440Myr gas t=440Myr 
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Fig. 20. Projected density of the model, at a time of 440 
Myr after the beginning of the collision: the stars Left , 
and the gas Right. The scale is in kpc. The black line is the 
trajectory of the companion starting from the southeast. 

The velocity fields do show the gross features of what is 
observed, to order of magnitude. In particular, they show 
reasonably regular rotation in spite of a quite perturbed 
galaxy "plane" . 



stars t=440Myr gas t=440Myr 
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Fig. 21. Same as Figure [20l for radial velocities. The 
color bar indicated at the right side goes from —210 to 
210kms-i. 

A feature that is well reproduced by the coUisional 
model is the sharp boundary in the southern edge of the 
gas disk (Figure [20|) . This sharp edge was one piece of 
evidence that had seemed to favor the ram pressure shock 
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scenario. An additional point of agreement is that the gas 
boundary appears in projection to be shifted to the north 
with respect to the stellar component boundary. This is 
because the gas disk is much thinner than the stellar disk, 
which is heated in the collision. The fact that the gas front 
was pushed "inside" the stellar disk was advanced as one 
of the main arguments to suggest a ram-pressure scenario 
by Shostak et al. p982p . In the present coUisional model, 
it is only a projection effect, arising from two components 
of different geometrical thickness. 

Another feature reproduced by the simulations is the 
" spoke" morphology connecting the two off-centered rings 
(see Sect. 4). This kind of feature is typical for head-on 
collisions, as is shown in the Cartwheel simulations (e.g. 
Hernquist & Weil, 1993; Horellou & Combes, 2001). These 
are visible in the models in the gas component only, simi- 
lar to the 8-micron dust-only image where they are clearly 
visible. Although the spokes in NGC 1961 are not at the 
same locations as those in our simulations, the similarity 
suggests that a head-on collision is capable of reproduc- 
ing, at least qualitatively, the main morphological features 
observed. 



6. Discussion 

The head-on collision model reproduces quite well the ring 
morphology and the sharp gas edge, which were difhcult 
to account for with ram pressure, given that hot IGM is 
almost absent in this small group of spiral galaxies. The 
minor merger scenario can also hold in the absence of a 
nearby interacting galaxy. 

One interesting outcome of the simulation is that it 
also gives a possible solution to the surprisingly large ra- 
dius of the galaxy with respect to its late (Sc) Hubble 
type. The outward propagating ring waves of the collision 
can increase the disk size to what appears to be abnormal 
proportions. 

One question that can still be asked is how such 
a galaxy in a relatively rich environment can retain 
such a late type, without a significant bulge. It may be 
that the galaxy is entering the group for the first time. 
Alternatively, even if it has already interacted and merged 
with some members of the group, they are all gas rich 
(Shostak et al. I1982p and could have replenished the 
disk. The in-plane accretion of a large number of satel- 
lite galaxies and intergalactic H I clouds with much lower 
masses can quietly build up the disk mass without dy- 
namically heating a spheroid (see Bournaud & Combes 
[2002)1 . thereby depressing NGC1961's bulge-to-disk ra- 
tio and driving its evolution towards later Hubble type. 
The most recent encounter has been sufficiently disrup- 
tive, however, that we would predict a reversion to an 
earlier type as the result of dynamical heating. The fact 
that we do not see more galaxies which are too massive for 
their Hubble types, suggests that this merger of moderate 
mass ratios must correspond to a quick transient phase in 
poor groups. 
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NGC 1961 is a low-luminosity AGN. From the central 
velocity dispersion, the mass of the central black hole can 
be estimated to be 3 x 10^ Mq (Krips et al. I2007p . From 
the X-ray luminosity of the nucleus and its deduced bolo- 
metric luminosity, its Eddington ratio is estimated to be 
-^boi/^EddlO^^ — 10"^. This means that the fuelling is not 
intense at this moment, and must be intermittent. The 
very different time-scales between dynamical phenomena 
at the extreme large and small scales is indeed an issue. 

NGC 1961 is however a good candidate for minor- 
merger induced fuelling of a low luminosity AGN; at least 
the minor merger made more potential fuel available in 
the inner kpc. This was proposed by several authors (see 
e.g., Taniguchi ll999p . and debated by others (e.g., Corbin 
I2000p . Hopkins & Hernquist (j2009p recently proposed to 
separate the fuelling of high-luminosity AGN (quasars), 
which would be fuelled by major mergers, from the low- 
luminosity ones, where secular evolution or minor merging 
is sufficient. NGC 1961 could be just an example of the 
transition between the two regimes. 

7. Conclusions 

Based on our observational and theoretical investigation 
of the asymmetries evident in NGC 1961 on both large and 
small scales, we have reached the following conclusions: 

1. The galaxy's large-scale optical asymmetries can be 
explained as the result of a head-on collision followed 
by a minor merger with mass ratio of ~ 1:4. Models 
that rely exclusively on ram pressure stripping cannot 
account for the level of disturbance observed in the 
galaxy's stellar light. 

2. Our preferred minor merger model can naturally ac- 
count for the excitation of m = 1 asymmetries in the 
galaxy's inner disk, and also for embedded off-centered 
rings in the CO maps. It can also explain the very 
sharp gas boundary in the southern part. Due to the 
geometrically thin gas plane, this sharp edge appears 
in projection shifted inside the stellar disk, a feature 
that had previously been thought to favor only a ram- 
pressure scenario. 

3. The head-on collision explains the ring wave propagat- 
ing outwards, but also the spokes with large pitch angle 
that connect the two rings. On the smallest scales, the 
relaxation time is sufficiently small compared to the 
merging timescale that the inner molecular disk has 
been able to become bar-unstable, as reflected in the 
m = 2 or nuclear bar morphology seen in the C0(2-l) 
map. The galaxy could have been barred before the 
merger, but the latter could have enhanced the bar 
strength, or create a new one, helping the gas inflow 
to fuel the AGN. 

4. Our collision scenario explains the exceptionally large 
size of NGC 1961 for its late type. The large disk warp 
accounts for an under-estimation of its true inclination 
in the plane of the sky, and thus an over-estimation 
of its dynamical mass. NGC 1961 may have built up 



a significant fraction of its current disk mass by the 
relatively quiescent accretion of intracluster gas. This 
scenario can account for the low bulge-to-disk ratio of 
this spiral galaxy in a relatively rich environment. The 
most recent and traumatic merger event, however, is 
likely to force NGC 1961 to revert to an earlier Hubble 
type. 

The simulations presented here only demonstrate the 
plausibility of the head-on collision, and arc not intended 
to provide an exact fit to the observations. In particular, 
the companion is only schematized as a rigid body, since 
there are insufficient observational constraints to define its 
physical parameters. We conclude however that NGC 1961 
could be a representative example of a minor-merger in- 
duced low-luminosity AGN. 
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